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Abstract 

o . 

, The Korteweg de Vries (KdV) equation with small dispersion is a model for the 

I formation and propagation of dispersive shock waves in one dimension. Disper- 

■ sive shock waves in KdV are characterized by the appearance of zones of rapid 

modulated oscillations in the solution of the Cauchy problem with smooth initial 
data. The modulation in time and space of the amplitudes, the frequencies and 
the wave-numbers of these oscillations and their interactions is approximately de- 
scribed by the g-phase Whitham equations. We study the initial value problem 
^ I for the Whitham equations for a one parameter family of monotone decreasing 

initial data. We obtain the bifurcation diagram of the number g of interacting 
oscillatory zones. 

1 Introduction 

The solution of the Cauchy problem for the KdV equation 

Ut + 6uUx + e'^Uxxx = 0, t,x eJR 
u{x, t = 0,e) = uo{x) 



:i.ii 



where e > is a small parameter, is characterized by the appearance of zones of rapid 
modulated oscillations. These modulated oscillations are called dispersive shock waves 
and they are described in terms of the Whitham equations. For e > 0, no matter 
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how small, the solution of ( |1 . 1| ) with smooth initial data exists and remains smooth 
for all if: > 0. For e = (|1 . 1| ) becomes the Cauchy problem for the Burgers equation 
Ut + 6uUx = 0. If the initial data uo{x) is decreasing somewhere the solution u{x,t) 
of the Burgers equation has always a point {xq, to) of gradient catastrophe where an 
infinite derivative develops. 



t>t„ 




Figure 1.1: The dashed line represents the formal solution of the Burgers equation 
after the time of gradient catastrophe t = Iq. The oscillations on the picture are close 
to a modulated periodic wave. 

After the time of gradient catastrophe the solution ^(a;, t,e) of (|1.1|) develops an ex- 
panding region filled with rapid modulated oscillations as shown in Fig. |1 . 1| . 

The idea and first example of the description of the dispersive shock waves were 
proposed by the physicists Gurevich and Pitaevski (GP) They considered the 
oscillatory zone to be approximately described by a modulated periodic wave: 



2a 



dn 



'6s 



Vt + xo) 



+ 7: 



(1.2) 



where dn{y,s) is the Jacobi elliptic function of modulus s G (0,1), xq is a suitable 



phase, the quantities a, s, 7 and V 



2a- 



3s 



+ 7 



depend on x and t. These 



quantities evolve according to the Whitham equations which guarantee the validity 
of the approximate description ( |1.2| ). For constant values of the parameters a, s and 7, 
u{x,t,e) is an exact periodic solution of KdV with amplitude a, wave number k and 
frequency ui given by the relations 



IT 



eK{s) 




a 
6s 



k_ 
V 



where K{s) is the elliptic integral of the first kind of modulus s. 
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Whitham introduced the Riemann invariants Mi > ^2 > W3 to write the equations for 
a, s and 7 in diagonal form. These quantities are expressed in terms of Ui > U2 > M3 
by the relations 

U2 — u?j 

a = U2 — U3 s= , 'y = U2 + U3 — ui. (1.3) 

Ml - M3 

The Whitham equations for the Wj, i = 1,2,3 read [0] 

dtUi{x,t) + Xi{ui,U2,U3)d^Ui{x,t) = 0, i = l,2,3, (1.4) 

where 

Ai(ui,U2,M3) = 2(ui + U2 + Ms) H , , (1-5) 

Ui + ao 



ao = -Ml + [ui - U3) 



K{s) 



and E{s) is the complete elliptic integral of the second kind. 

The variable U2 can vary from M3 to Ui. The oscillatory region is bounded on one side 
by the point x~{t) where U2{x,t) = U2,{x,t), and on the other side by the point x~^{t) 
where U2{x,t) = ui{x,t) (see Fig. and Fig. |1.2| ). Outside the region {x~ (t) , x~^ (t)) 
the solution u{x, t, e) of ( |1.1D is well approximated by the solution u{x, t) of the Burgers 
equation. 

In the theory of singularity a function behaves quadratically in the neighborhood of its 
generic singular point |^ . In a similar way, in the theory of dispersive shock waves the 
generic analytic monotone decreasing initial data x = f {u\t=o) behaves (up to shifts 
and rescalings) like x = —u^ in the neighborhood of its generic breaking point. Thus 
the behavior of the generic analytic initial data near a point of gradient catastrophe is 

X - Xo ^ -C3(m - Uof - Ci{u - UqY Ck{u - Mq)^ " (1-6) 

GP kept only the first term of the above Taylor series namely they considered the 
initial data x — xq = —{u — uqY (up to rescaling) and solved the equations (|1.4| ) 
numerically for such initial data. Potemin obtained the analytic solution of the 
equations ( |1.4D for the same initial data and showed that x^{t) = xq — 12\/3 (t — tg)^ 



and x+(t) = xo + 4/3 v^VS (t - to)^ (see Fig. [Ol) . 



The higher order terms of the Taylor expansion (|1.6|) are negligible only for t ^ to + 

one starts considering the higher order terms of the expansion 
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, u(x,ti) 



g=o 



"^Ul (x,ti) 



x"(t,) 



a) 



t=t, 




g=o 



~ u(x,ti) 



x-(t) 



g(x,t)=0 



x-(t,) 



:(x,t)=l 




b) 



:(x,t)=0 



x*(t,) 



Figure 1.2: On picture a) the dashed hue represents the formal solution of the Burgers 
equation, the continuous line represents the solution of the Whitham equations. The 
solution {ui{x,t) , M2(x, t), U3(x, t)) of the Whitham equations and the position of the 
boundaries x~{t), x~^{t) of the oscillatory zone are to be determined from the condi- 
tion u{x~{t),t) = Ui{x~{t),t), Ux{x~{t),t) = Uij:{x~{t),t), U{x^{t),t) = U3{x^{t),t), 

Ux{x~^ (t),t) = U3^(x+(t),t) where u{x,t) is the solution of the Burgers equation. Pic- 
ture b) represents the functions x~{t) and x~^{t) on the x — t plane. 



( |1.6|) , we come to a situation in which two or more oscillatory wave trains come in 
interaction or in the case in which the solution of (|1.4|) comes itself to a point of 
gradient catastrophe. We will see that such a phenomena necessarily occur in the 
simplest non trivial deformation of the cubic law which turns out to be the fifth order 
polynomial: 



X 



-C3(m - Uof - Ci{u - UoY - C5(u - Uq)^ , C3 > 0, C5 > 



(1.7) 



Through the shift (x x+6tuo, u — > u+Uq) the above initial data becomes x ~ —c^u^— 
c^u^ — c^u^ . The parameter space (03, C4, C5) can be reduced exploiting the invariance of 



the KdV equation under the groups of transformations 



X 



PL 



u 



ku) 



k and {x ^ ax, t — > at), a 7^ 0. These transformations change, however, the 
value of the small parameter e. Taking k = ^/cs/c^ and a = C3, the initial data (|1.7] ) 
can be reduced to the form 



X = —u 
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where the dimensionless parameter c is chosen in the form 



V5 no 



The monotonicity condition requires < 15/4. From the above considerations the 
polynomial ( p. .81) represents the generic one-parameter deformation of a monotone de- 
creasing analytic initial data with cubic inflection point. 

For cubic initial data the solution u{x, t, e) of (|1.1| ) is approximately given by equation 



(|1.2|) inside the oscillatory zone and by the solution of the Burgers equation outside this 
zone. For general initial data one have to use higher genera analogues of the Whitham 
equations (|1.4|) to approximately described the solution u{x,t,e) of (|1.1|) 0,0 • Lax 
and Levermore [Q, and Venakides ^ described, for certain particular classes of initial 
data, the dispersive shock waves in the frame of the zero-dispersion asymptotics for 
the solution of the inverse scattering problem of KdV. According to their results, to 
the solution u{x,t,e) as e ^ it corresponds a decomposition of the {x,t) plane into 
a number of domains Dg, g = 0,1, ... . In the domain Dg the principal term of the 
asymptotics is given by the g-phase solution of the KdV equation 

u{x, t, e) = $ f^^lM, . . . , u,{x,t),..., u,g^,ix, t)] , (1.9) 



where the functions Sj{x,t) satisfy the equations JTO 



dS- dS- 

= kj{u{x,t)) , = ujj{u{x,t)), j = l,...,g, (1.10) 

and the formula 

u{x, t) = $ {kix + cjit + (/>!,... , kgX + ujgt + (t)g; Ml, ... , U2g+l) 

for constant values of the parameters Mi, ... , U2g+i, kj = kj{u) and Uj = ojjiu) and for 
arbitrary 0j, j = 1, . . . ,g, gives the family of the so-called g-phase exact solutions of 
KdV for e = 1 [Q. The wave parameters in ( p..lO| ) depend on the functions ui{x,t) > 



> U2g+i {x, t) which satisfy the system of equations 

+ Xi{ui,U2,. . . ,U2g+i)-^ = , i = l,...,2g + 1 , g>0. (1.11) 



For a given g the system { \L . 1 1|) is called (/-phase Whitham equations. The Whitham 



equations is the collection of all systems (|1 . 1 1| ) for g > 0. The zero-phase Whitham 
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equation coincides with the Burgers equation. For g > the speeds Xi{ui,U2, ■ ■ ■ , M2g+i), 
i = 1, 2, . . . , 2(yf + 1, depend on mi, . . . , U2g+i, through complete hyperelhptic integrals of 
genus g. For this reason the (/-phase solution is also called solution of genus g. When 
the (7-phase solution comes to a point of gradient catastrophe, it should be continued 
by the {g + l)-phase solution. The main problem is to glue together the solutions of 
the Whitham equations for different g in order to produce a C^-smooth curve in the 
X — u plane evolving smoothly with t . This problem is referred to as initial value 
problem of the Whitham equations. The genus g{x,t) is a piecewise constant on the 
(x, t) plane (see Fig. |1.2| ). For generic initial data it is not not whether the genus g{x, t) 
is a bounded function. 

The algebraic geometric description of these equations for (7 > 1, was first derived in 
|p[. Dubrovin and Novikov developed a geometric Hamiltonian theory for the Whitham 



equations [|T2]- Based on this theory, Tsarev |jTT[ was able to prove that, for each 
g, equations ( |1.11D can be solved by a generalized method of characteristic. This 
method was put into an algebro-geometric setting by Krichever [|12| . Recently Dubrovin 
proposed a variational principle for the Whitham equations. In the minimizer of 
a functional defined on a certain infinite-dimensional space formally solves the initial 
value problem of the Whitham equations for each point of the (x, t) plane. 
We use this variational principle to arrive at the main result of this paper. This is a 
complete description of the bifurcation diagram in the x — t plane of the one-parameter 
family of initial value problems for ( |1 . 1 1| ) . For the purpose we use the result in 
where it is shown that for polynomial initial data of degree 2N + 1 the solution of 
the Whitham equations has genus at most equal to A^. In our case of fifth degree 
polynomial initial data, the solution of the initial value problem turns out to be glued 
from solutions of the Whitham equations of genera g = g = 1 and g = 2. The 
various topological type of bifurcation diagrams are described in terms of the number 
of points of gradient catastrophe that appear in the solution of the Whitham equation 
and other particular points which we call double-leading edge, double trailing edge and 
leading-trailing edge. 

This paper is organized as follows. On Sec. 2 we give some basic results about the 
theory of the Whitham equations and we describe Dubrovin's variational principle for 
the Whitham equations. In Sec. 3 we describe the phase transition boundaries for 
general initial data. In Sec. 4 we study the function g = g{x,t) for the one parameter 
family of initial data ( |1.8|) . 
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2 Preliminaries on the theory of the Whitham equa- 
tions 

The (7-phase Whitham equations ( p, . 1 1| ) are constructed in the following way. 
On the Riemann surface 

r 2g+l ■\ 

Sg := < (P, /u(r)), /U^ = J]^ (r - Mi) > , Mi > M2 > ■ ■ ■ > ^23+1, Ui G IR, (2.1) 

we define the abelian differentials of the second kind Ck, k > with a pole at infinity 
of order 2/c + 2 



(2.2) 



where u = (wi, . . . , M23+1) and the constants ai = ai{u), i = 0,1, . . . , k + g — 1 , are 
uniquely determined by the conditions: 



al{r) = 0, j = l,2,...,g, al{r) 



rU2j 

We use the notation 



2 



+ 0(r 2) 



dpg{r) = (T^(r), dqg{r) = 12af (r) . 



dr for large | r | (2.3) 



(2.4) 



In literature the differentials dpg{r) and dqg{r) are called quasi momentum and quasi- 
energy respectively [|r^. The speeds Xi{u) of the g-phase Whitham equations ( |1.11| ) 
are given by the ratio H: 



Xi{u) 



dqg{r) 



dpg{r) 



i = l,2...2g + l. 



(2.5) 



dr 12r — 6u 

In the case g = 0, we have dpo[r) = —= and dqo[r) = — dr respectively, 

y/r — u y/r — u 

so that the zero-phase Whitham equation (|1.11|) coincides with the Burgers equation 
Ut + 6uUx = . We consider monotone decreasing polynomial initial data of the form 



X = f{u) = Co + CiU-\ h C2NU^^ + C2N+1U 



,2N+l 



(2.6) 
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For such initial data the solution of the Burgers equation is obtained by the method 
of characteristic and is given by the expression 



X 



6tu + f{u) 



(2.7) 



This solution is globally well defined only for < t < tc where the time tc = 
I minug]R[— /'(u)] is the time of gradient catastrophe of the solution. The breaking 
is cause by an inflection point in the initial data. For t > tc the solution of the 
Whitham equations is obtained gluing together C^-smoothly solutions of different gen- 
era. The (/-phase solution ui > U2 > ■ ■ ■ > ii2g+i which is attached smoothly to the 
{g + l)-phase solution is obtained by a generalization of the method of characteristic 
and reads 



X = Xi{u) t + Wi{u), Wi{u) — ^ 



dpg{r) 



l,...,2g + l 



(2. 



where Xi{u) has been defined in ( |2.5|) and the differential dsg{r) is given by the expres- 
sion 



2N+1 



fc=0 



(2.9) 



where cr^ has been defined in 



and the c^'s have been defined in ( p.6|) . 
The solution (|2.8| ) can be written in the equivalent algebro-geometric form. Consider 
the differential 



fi,(r) 



Xdpg 



tdqg{r) - dsg 



where dpg, dqg and dsg have been defined in (|2.4| ) and (| 



(/-phase solution is given by the relation |I2| , |T7| 



ilgir) 







1,2. ..,2(7 + 1. 



(2.10) 

respectively. Then the 

f2.in 



The solution ui > U2 > ■ ■ ■ > M2g+i of the (7-phase Whitham equations ( |1.11D is 
implicitly defined as a function of x and t by the equations ( |2.8|) or ( |2.11| ). The 
solution is uniquely defined only for those x and t such that the functions Ui{x,t) are 
real and d^Ui^x, t), i = 1, . . . ,2g + 1, are not vanishing. 
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2.1 Variational principle for the Whitham equations 



The solution of the Whitham equations, for a given initial data, can be written as the 
minimizer of a functional defined on a certain infinite-dimensional space [l^. Let us 
first consider the zero-phase equation. The characteristic equation x = 6tu + f{u), 
where f{u) is the polynomial (|2.6| ), can be consider as the minimum of the function 



xu 



3tu^ - F{u) 



(2.12) 



where x, t and c = (cq, Ci, . . . , C2Ar+i) are parameters and F'{u) = f{u). The minimum 
is unique only for t < tc- In the function of type ( |2.12| ) is extended to a functional 
onto the moduli space of all hyperelliptic Riemann surfaces with real branch points 
in such a way that the absolute minimum of this functional gives the solution of the 
initial value problem of the Whitham equations ( |1 . 1 1|) for the initial data ( |2.6| ). 

First we define the restriction of this functional on the hyperelliptic Riemann sur- 
faces of genus g with real branch points Mi > M2 > ■ ■ ■ > M2g+i- The restriction is a 
function depending on the branch points. This function is build as follows. Consider 
the asymptotic expansion of the quasi- momentum dpg (see e.g. [||) 



dpg{r) 



^ V fc=0 



dr . 



(2.13) 



The coefficients = /^('Ui, ^2, • • • , ^23+1), are the so called KdV integrals and are 
smooth functions of the parameter Mi > ^2 > ■ • • > M2g+i- 



Theorem 2.1 JT^/ On the Riemann surface Sg consider the function 

2N+1 



Gl,t,c]i^) = -^Io{u) + 3th{u) - J2 



k\ 



k=0 



2^'(2A; - 1) 



-Xkhiu) 



(2.14) 



depending on the variable ui, U2, ■ ■ ■ , ^29+1 one? on the parameters x, t and c. Then the 
equations 



_d_ 

dui 



Gl^Au) = Q, z = l,2... ,2^7 + 1 



(2.15) 



are equivalent to the equations ^g{r) \r=u,= , i = 1, . . . ,2g + 1, where ^g{r) has 
been defined in 
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Remark 2.2 The function G|^^^j('u) can also be written in the form 

Gf^.Au) = Res[dp,iz)iJ^iz) - 2^zx)i T{z) = f ^^^±Md^, (2.16) 

L ' ' 1 z=oo - ^ 

and Res2=oo is the residue evaluate at infinity. 

To extend the function ( p.l4|) defined on the hyperelhptic surfaces Sg>o to a functional 
on the infinite dimensional space M of all hyperelhptic Riemann surfaces Sg>Q and 
their degeneration, we follow |jl3| . 

Construct the space M inductively starting from Mq = IR. We denote u the coordinate 
in Mq. Define now 

Mg = M° U Vi U T,_i , 

where 

M° = {{Ui,U2, . . .,U2g+l) e IR^^+^I Ui>U2> ■■■> U2g+l} , 



Each space -f'^-i is attached to the component of the boundary of M° where 

U2j-i-U2j^0, j = l,2,...,g; 
each space is attached to the component of the boundary of M° where 

U2j - U2j+i , j = l,2...,g. 

For reasons that will become clear later, we call leading edge each boundary component 
j = ^, ■ ■ ■ , 9 and trailing edge each boundary component Tg_-^ j = 1, . . . , g. 
A smooth functional / on M is a sequence of function fg defined on every Mg and 
having certain good analytic properties on the boundaries of the space Mg (for details 
cfr |TB[). These properties are illustrated in the following example. 
Example. Let us consider the Riemann surface of genus g + 1 



25+1 

n 



fi^ = {r — V — y/5){r — v + \/5) Y\_(^ ~ ^i)' ^ ^ IR-, ui > U2 > ■ ■ ■ > U2g+i, 
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where v ^ Uj, j = 1, . . .2g + 1 and < 5 ^ 1. We denote with dpg-^.i{r,S), the 
quasi momentum defined on this surface and with G^^^^-^{u,v,6) the function equiv- 
alent to (|2.14|) defined on this surface and depending on the variables v — ^/5, v + 
Vd, Ml, ... , U2g+i- We study the behavior of dpg+i{r, 6) in the limit 5 — >• 0. 
First we consider a trailing edge of the boundary of the space M^+i, namely the case 
in which v G {u2j, M2j-i), 1 < j < + 1, M2g+2 = —oc. We have that [|T8] 

rfpg+i(r,5) = dpgir) + ypg{v)^uj,{r) + 0i6'') (2.17) 

dp (t'^ 

where dpJv) = — — \r=v and 0(6'^) /S"^ is a meromorphic differential with zero residues. 
dr 

The differential ujy{r) is a normalized Abelian differential of the third kind with poles 
at the points P^{v) = {v, ±fi{v)) with residue ±1 respectively, namely 

u;^ir) = ^^-j2"-^B,iv), /i^) = - • (2.18) 

fi{r}r-v ^ /i(r) ^J. 

The constants Bk{v) are uniquely determined by the normalization conditions 

uj^{r) = 0, k = l,...,g. 

When V G {u2j+i,U2j), j = 0,...,g, and uq = +oo, we are considering one of the 
leading edge of the space M^+i. 



In this case the expansion of the differential dpg+i{r, 6) for 5 — > reads |T8 

dpg+,ir, 6) = dpg{r) + puj,{r) / dpg{0 + O^p") (2.19) 

Jp-{v) 

where ujy{r) has been defined in (|2.18| ) and p = The correction in the right hand 

side of (p.l7|) contains also exponentially small terms like exp(— 1/p). 
Combining equations ( |2.11| ), ( [^.16D and ( ^.17| ) we obtain the expressions of the function 
G^^^^-^{u,v,6) near the trailing edges of the space M^+i, namely 

Gg,j(M,M,5) = Gf.,,,](^) + yPgiv)Qgiv) + Oi6'), (2.20) 
where the differential ^g{r) has been defined in ( p.lOj ) and here and below ^g{v) = 



dr 
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The behavior of the function G^^^^^{u,v,S) near the leading edges of the space M^+i, 
is obtained from ( pTTD , {^A§ and (^J^): 



GlS^u, V, 5) ^ Gf,,,„(«) + dp,{r) r n,{r). {2.21] 

^ Jp-{v) Jp-{v) 



The function ^ ^ j (wi, 1*2, ... , "^29+1) defined on M° can be extended to a smooth func- 
tional on the space M. In the extension is built proving that the /^(mi, . . . , M2g+i), 
> 0, can be extended to smooth functionals on M. We state the following theorem 



Theorem 2.3 /[7^ T/ie functional 



k\ 



k=0 



is a C°° smooth functional on M . Its absolute minimum is a -smooth multi-valued 
function of x depending -smoothly on the parameters t,ci, . . . , C2n+i- 
If the absolute minimum {ui{x,t), . . . ,U2g+i{x,t)) belongs to M° for certain values of 
the parameters, then it satisfies the g-phase Whitham equations. 



3 Study of the function g = g{x, t) 

In this section we study the behaviour of the genus g = g{x,t) of the solution of the 
Whitham equations (|1 . 1 1|) for the one-parameter family of initial data (|1.8|) . In |jl4| 



it was shown that for polynomial initial data of degree 2N + 1 the solution of the 
Whitham equations has genus at most equal to A^. Hence in the case of initial data 
(|1.8| ) the solution will have genus at most equal to 2. For this reason we need only to 
describe the locus of points of the x — t plane where the genus g{x,t) increases from 
zero to one or two and from one to two. To this end we study the behaviour of the 
functional ( p.22|) on the phase transition boundaries of the space Mq, Mi and M2. The 



restriction of the functional (p.22|) on Mq has the form (|2.12| ). Thus the solution of the 



Burgers equation (i.e. the Whitham equation for g = 0) 

x = 6tu + f{u) (3.1) 



minimizes the functional ( p. 221) for all x and t < t^ where tc is the time of gradient 



catastrophe. Let us find for which x and t the solution of (|3.1| ) is also minimal along 
directions transversal to Mq C M. For the purpose we consider the embedding of Mq 
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as the component of the boundaries of Mi (To, Lq) and we study the behaviour of the 
functional (|2.22|) near the traihng edge and leading edge of the space Mi. 
First we need the following theorem. 

Theorem 3.1 On the solution ui{x,t) > U2{x,t) > ■ ■ ■ > U2g+i of the g-phase 
Whitham equations, the differential ^g{r) defined in (^.ICi) reads 



2g+l 



n,(r) 



-fi{r)[dr-^g{r, u) + J2 9u,'^9ir, u)]dr 



(3.2) 



i=l 



where the function \E'g(r, u) is given by the expression 

J-'{z)dz 



^g{r, u) = — Res 



_li{z){z - r)_ 
and Res^=oo is the residue evaluate at infinity 



, ^(z) 



6te + /(0 
Vz-^ 



(3.3) 



Using ( 2.2C ), (|3.2| ) and ( p.3|) we can simplify the functional (|2.22 ) near the trailing 
edge of the space Mi to the form 

GUtd^,v,6) = Gl,^,^iu) + ^(-6t - dMu.v) - d^qiu,v)) + 0(5') , (3.4) 



where v < u, G? ^ ^, (m) has been defined in ( p.l2|) and 



q{u, v) = — Res 



{z — v)\/ z — u 



1 /(l±!Zi^ + 1^^) 

^ dm. 



4^2 7-1 



m 



(3.5) 



Analogously on the leading edge of the space Mi we obtain 

G\,^,^4u,v,b) = GjL,,,,](w) + 16p{v - u)\-2t - duq{u,v)) + O(p') . (3.6) 

If, for fixed (m, t), the (5-correction of ( p. 41) is positive for every v < u, then the minimizer 
belongs to Mq. If it is negative for some values of v < u, then the minimizer belongs 
to the inner part of Mi. The points v < u belong to the trailing edge of the space Mi 
if the triple (t, m, v) is a zero and a minimum with respect to v of the 5 correction of 
(|3.4|) . From these considerations we obtain the following lemma. 
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Lemma 3.2 The points v and u, v < u, belong to the boundary Tq of the space Mi if 
(t, u, v) satisfy the system 



Qt + dvq{u, v) + duq{u, v) = 

d^{dyq{u, v) + duQiu, v)) = (3.7) 
{d^yid^qiu, v) + duq{u, v)) <Q. 

The curve x = x~{t) on the x — t plane where the genus increases from zero to one 
is determined solving the above system together with the zero-phase solution x = 
6tu + f{u). We will call the curve x~{t) trailing edge. 

Analogous consideration can be done for the leading edge. 

Lemma 3.3 The points v and u, v > u, belong to the boundary Lq iff {t,u,v) satisfy 
the system 

-2t + duq = 

OAq = (3.8) 
{d,yduq<0. 

The curve x = x'^{t) on the x — t plane where the genus increases from zero to one 
is determined solving the above system together with the zero-phase solution x = 
6tu + f{u). We will call the curve x^{t) leading edge. 



Remark 3.4 Both systems 1\ ) and 8^ ) in the limit v ^ u become 

m + f'{u) = 

fiu) = (3.9) 
riu) < 



where f{u) is the initial data ( ^.dj ). If the above system admits a real solution {u,t) 
such that 

2t + d^q<0 Wvj^u 

6t + dyq + duq<0 yvy^u, ^ ' 

thenu belongs to the boundary T^PiLo of the space Mi. The corresponding point {x,t,u) 
is a point of gradient catastrophe of the solution of the Burgers equations. 
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3.1 Boundary between the one-phase solution and the two- 
phase solution 

In the following, we write explicitly the equations determining the phase transition 
boundary between the one-phase solution and the two phase solution. 
Let be dpi{r) the quasi- momentum restricted to the inner part of Mi, 

r"2 rdr 

(r) - ^ + "0 a - ' V(--"i)(--»^)(--"3) , . 

api[r) - ==ar, ao - — 7^75 l-^-iiJ 

^/{r - ui){r - U2){r - U3) J , 

From (p.l6| ) the restriction on the inner part of Mi of the functional G[x,t,c] reads 

<^Mci(«i)«2,M3) = Res[dpi{z){J^{z) - 2^/zx)]. (3.12) 



Proposition 3.5 The critical points of ( \3.12{ ) on the space of elliptic curves are given 
by the equations for ui > U2 > U3, 



X = Xit + Wi , z = 1, 2, 3 



where A, = Xi{ui,U2,us) has been defined in and 



n jj^i ~ 



Wi[Ui,U2,U3) 



aQ + u. 
a{ui, U2, Us) = - Res^=oo[ 



-du^a{ui,U2,U3) 



Jo 



d^ ) dz 



Proof: Using the formula |13 



dui ° 2 2nj^i(Mi-MjO 



z= 1,2,3, 



(3.13) 



(3.14) 



(3.15) 
□ 



equations ( p.l3| ) and (|3.14|) are recovered straightforward. 
The trailing edge of the phase transition boundary between the one-phase solution and 
the two-phase solution is determined from ( |2.20|) , ( pl2|) and ( p.3|) . When g = \ (|3.3| ) 
reads 



^'i(r, u) = — Res 



T{z)dz 



{z - r) ^J{z- ux){z - U2){z - U3) 



(3.16) 
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so that ( p.20| ) simplifies to the form 



GUt,c]i^, V, S) ^ Gl^,^,]{u) + -{v + ao){-d,^i{v, u) -J^ du^'^ii^, ^)) (3-17) 



i=l 



where Gl 



[x,t,c] 



has been defined in ( p.l2| ) and ao has been defined in ( p. 11 ). The 
point V G (— 00,^3) or f G {u2,Ui). In this case v + ao ^ 0- From the 6 correction of 
(|3.17|) we obtain the equations determining the traihng edge, namely 



[ i^,mv + ao){^,<^I,{v,u) + J2l,^^^<^I^{v,u))<0 



(3.18) 



where u = {ui{x , t) , U2{x , t) , u^i^x , t)) is determined from ( |3.13| ). A similar system can 
be obtained for the leading edge using ( |3.6|) . 



In the limit v — >• u;, 1 < Z < 3, system ( p. 181 ) reads 
dut{du^<y + du2(T + du.^a) = 

(duifidu^a + du^cr + du^a) = 



(3.19) 



where a = cx^ui, U2, Ms) has been defined in ( p.l4|) . If the above system admit a solution 
x*,t*, then the one-phase solution ui{x*,t*) > U2{x*,t*) > U3{x*,t*) determined by 

d 

(13.131) has on the ui branch a point of gradient catastrophe, namely —ui(x, t) 
ox 



x — x 

t=t* 



00 and 



dx'^ 



ui{x,t) 



X — X^ 

t=t* 



00. 



Theorem 3.6 // the initial data ^2.(^ ) satisfies the condition f"'{u) < 0, then the 
solution of the one-phase Whitham equations has no point of gradient catastrophe. 



Proof: Using the following expression for cr(Mi, U2, 



0-(Ui,M2,M3) 



2^/271 J -I J -I 



¥-3) 



dtds, (3.20) 
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it can be easily checked that the second equation of ( p.l9| ) cannot be satisfied for any 
real Ui,U2 and M3. □ 



Remark 3.7 Using a different approach Tian has proved j^T^that for smooth monotone 
decreasing initial data satisfying the condition f"'{u) < 0, the solution of the one phase 
Whitham equations exists for all t > 0. 



4 Bifurcation diagrams 



We study the bifurcation diagram of the solution of the Whitham equations for the 
one parameter family of initial data 



fciu) 



For such initial data, the functional (|2.22| ) reads 
G 



\x,t,c\ 



Xin + 3t/i ch 

u ^ 20 70 



252 



(4.21) 



(4.22) 



The restriction of this functional on Mq (that is on the curve /z^ 



the form 



u 



G'[^,t,c](w) =xu-?>tu + — + C— + 



u, u eJR) has 



(4.23) 



thus the minimizer given by 



X = Qut — u — cu — 



(4.24) 



solves the Burgers equation until the time of gradient catastrophe to = 0. At later 
times the minimizer of G\x^t,c] may belong to Mq, Mi or M2. 

The Burgers equation has another point of gradient catastrophe if ( p.9|) and ( 3.10 ) with 
the initial data ( [4.21|) are satisfied, that is 

- 6t + 3^2 + Acu^ + = 
6m + I2cu^ + 20^3 = 
6 + 24cM + 60^2 > 

6t + d^qc + duqc <0 ^ v ^ u , fGlR, 
^ 2t + duqc <0 Vt;^u, wGlR, 



(4.25) 
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where 



qc{u, v)=- 1/35 {5u^ + 6u^v + 8uv'^ + 16v^) - c/315 (35 + 
AOu^v + 6Auv^ + 48 t;^ + 128 v^) - 1 /693 (63 + 
70 M% + 80 + 96u'^v^ + 128uv^ + 256 v^) . 



The solutions of (|4.25| ) are obviously Uq = 0, to = and 



Ml 



-— c + 

10 ^ 
10 



for c > 



f) forc<0 



^1 = -Toic- 

ti = U^uj + Acuf + 5uj) > 



with the constraints 



4^<c<-JV5(13 + 5v^), y^(75 + 21v^)<c<^ 



which are obtained from the last two inequalities of ( 4.25 ). We put 

1^1 = + t)VY) ^ -1.90863 

Z/4 



-^^5(13 + 5^7) : 

— (75 + 21v^) ^ 1.88494 
44 ^ 



(4.26) 



(4.27) 



(4.28) 

(4.29) 
(4.30) 



because these numbers occur frequently in the following. The xi coordinate relative to 
[ti, Ui) is recovered from the equation (|4.24| ). For c positive Xi > and for c negative 
xi < 0. 

Here and below all the numerical results are obtained using Mathematica 3.0 for Solaris 
Copyright 1988-97 Wolfram Research, Inc. 



4.1 Trailing edges 

The equations determining the trailing edge of the space Mi are given by system (|3.7|) , 
namely 

6t + d^Qc + duQc = 

d,id,q, + duqc) =0 (4.31) 
{d^Yid^qc + duqc) < , 

where qc{u,v) has been defined in ( [4.26| ). The above system can be easily solved 
numerically getting u = u{t) > v = v{t). Substituting the value of u = u{t) in the 
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zero phase solution ( [4.24| ) we get, in the x — t plane, the curve x^{t) which describes 
the phase transition from the g = solution to the g = 1 solution. For some values of 
c and t system (|4.31|) has two different solutions ui{t) > vi{t) and U2(t) > V2{t) and 
correspondingly we have two curves xj"(t) and X2 (t) in the x — t plane. This occurs 
for c in the intervals 



15 



< C < Z/i, Z/4 < C < 



where and z/4 have been defined in ( ^.29| ) and ( [4. 301 ) respectively. 

A particular solution of system (|4.31| ) is the one which describes the double-trailing 
edge (cfr. below). System ( 4.31 ) determines the zeros which are also a maxima with 
respect to v of the polynomial 6t + d^qc + duQc- Since d^qc + duQc is a fourth degree 
polynomial in the variable v with negative leading coefficients, it cannot have more 
than two maxima. When system ( [4.31| ) admits such two maxima (t, u, v) and (t, u, w) 
with i > 0, V < u and w < it, we are in a situation of a double trailing edge (see 
Figure JTT). The degenerate Riemann surface describing this situation is given by the 
equation /i^ = (r — u){r — v)'^{r — w)'^ and represents the boundary T/ fl of the 
space M2. The curve x~{t) describing the trailing edge in the x — t plane looses the 
C^-smoothness at the point t = t. Such two solutions exists for c in the intervals 



1 D. 




Di 


D„ 


u(xj~ — 




u,(x) 






u,(x) 






U3(x)=U2(x)=W^ 






















X 










U5(X)=U4(X)=V/' 






















-^u(x) 



Figure 4.1: Double trailing edge. 



15 V15 
< C < Z^i , z/3 < c < 
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where ui has been defined in ( |4.29| ) and 



1^3 




(^(27 - 7^21)^ + (27 + 7^21)^) ~ 1.78167 . (4.32) 



For c = the two solutions (t, u, v) and (i, u, w) become coincident, namely v = w. 
For c = i^i the solution (t, it, w) satisfies w = u. The point {v, w = u) belongs to 
the component L\ fl Tl fl of the boundary of M2 and the corresponding degenerate 
Riemann surface is /x^ = (r — uY[r — vY- The solution of the Burgers equation has in 
(x, i,w = u) a point of gradient catastrophe in correspondence of the trailing edge. 
For c in the interval ui < c < there exists just one real solution v{t,c) < u{t,c) of 
system ( [4.3 1| ) for all t > 0. Substituting u{c,t) in the zero-phase equation we get a 
curve x~{t) in the x — t plane which is smooth for all t > 0. 

4.2 Leading edges 

The leading edges of the space Mi are determined from system ( p.8| ) namely 

2t + duQciu^v) = 

dyduQciu^v) = (4.33) 
{dyfduqc{u,v) < 



where qc{u,v) has been defined in ( |4.26| ). For some values of t > and for c in the 
intervals 



15 V15 
^<c<z/i, z/4<c<— ^ 

system ( [4.33| ) has two different solutions ui{t) < vi(t) and U2{t) < V2{t). Correspond- 
ingly we have two leading edges Xi{t) and x^(t) in the x — t plane which describe the 
phase transition from the g = solution to the g = 1 solution. 

We also consider the situation of a double leading edge as shown in Fig. |4.2| . That is we 
study for which values of c system ( |4.33| ) has two real solutions {t, u, w) and (t, £t, v) 
satisfying the constraints v > u, w > u and t > 0. These solutions belong to the 
boundary component L\ fl of the space M2 which is described by the degenerate 
Riemann surface fi"^ = {r — u){r — {;)^(r — w)^. This double solution occurs only for c 
in the intervals 

/15 /l5 
20 



where z/4 has been defined in ( [4 .301 ) and 



/5 35 TT-^ ^11 , , 

For c = 1^4, these two solutions become coincident, namely v = w. 

For c = z/2 the solution (t, u, w) satisfies the relation w = u. The corresponding 
point (x, t, w = £t) is a point of gradient catastrophe of the Burgers equation and it 
occurs at the leading edge. 

When c G (z/2, 1/4), there exists just one real solution f (t, c) > u{t,c) of system ( [4.31| ) 
for all t > 0. In this case there exists just a single leading edge for all t > 0. 




Figure 4.2: Double leading edge. Figure 4.3: Leading- Trailing edge. 



4.3 Leading-trailing edge 

We call leading-trailing edge the situation in which a leading edge and a trailing edge 
have the same [x, t) coordinates as shown in Fig. |4.3| . In other words a leading trailing 
edge is determined by a solution (t, u, w) of system ( [4. 331 ) and (t, -u, v) of system ( [4.31| ) 
where w > v. We check numerically that such two solutions exist for c in the intervals 
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4.4 Point of gradient catastrophe of the one-phase solution 



The solution of the one-phase Whitham equation with the initial data ( [4.21| ) has a 
point of gradient catastrophe if the correspondent system (|3.19|) has a solution. From 
Theorem ^]6| it follows that for < | , the solution of the one-phase Whitham equations 
with initial data ( 4.21 ) has no point of gradient catastrophe. 

We solve numerically systems ( |3.19|) a nd ( p.l3|) fo r the initial data ( |4.21| ) restricting 
c in the intervals -v^/2 < c < -y/5/2 and ^jijl <c< ^15/2. 

We find that there exists a point of gradient catastrophe for the solution of the 
one-phase Whitham equations on the -Ui-branch for c G [z/i, z/2), where Ui ^ —1.90863 
and i>2 — —1.85585 have been defined in ( [4. 291 ) and ( [4.34| ) respectively, (see figure |4.4] ). 

There is a point of gradient catastrophe in the solution of the one-phase Whitham 
equations on the M2-branch for c in the intervals [—^^, 1^2) and (2/3, ^^), see figure ^A . 

There is a point of gradient catastrophe in the solution of the one phase Whitham 
equations on the M3-branch for c G (z/3, z/4], where z/3 ^ 1.78167 and z/4 ~ 1.88494 have 
been defined in ([4. 291) and (|4.3CI|) respectively, see figure O. 



Remark 4.1 A point of gradient catastrophe appears in the solution of the one-phase 
Whitham equations whenever the solution is changing genus. If the point of gradient 
catastrophe appears on the Ui or u^-hranch, the corresponding solution will increase 
genus by one near this point. 

If the point of gradient catastrophe appears in the one-phase solution on the U2-branch 
it means that the two-phase solution has just disappeared. We have checked numerically 
this fact solving system l \3.18{ ) in a neighborhood of each point of gradient catastrophe. 
System ( ^.ISj ) determines the phase transition boundary between the one-phase solution 
and the two-phase solution. When a point of gradient catastrophe {x*,t*) appears on 
the Ui or u^-branch, then system ( ^.ISj ) admits a solution for t > t* . Namely a two- 
phase oscillatory zone is developing. On the contrary if a point of gradient catastrophe 
{x*,t*) appears on the U2 branch, then system (\3.1d[ ) does not have admissible solutions 
for t > t*. 



4.5 Bifurcation diagrams in the x — t plane 

We draw in the (x, t) plane the various topological types of bifurcation diagrams of the 
solution of the Whitham equations with initial data x = —v? — cu^ — , < ^• 
We draw with a solid line the points of the (x, t) plane where the genus increases from 
zero to one and with a dashed line the points where the genus increases from one to 
two. We have several cases. 
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1) i>'4 < c < Y ^; there exists a second breakpoint for the zero-phase solution in xi > 0, 

^1 > 0; there exists a point of gradient catastrophe on the it2-branch of the one-phase 
sohition for t > ti, there are a double trailing edge, a leading-trailing edge and a double 
leading edge hence the bifurcation diagram of the genus g{x, t) is 



point of gradient catastrophe of the one-phase solution 








2a) c = z/4; the Burgers equation has a second point of gradient catastrophe in corre- 
spondence of the leading edge, there is a point of gradient catastrophe on the M2-branch 
of the one-phase solution and there is a double trailing edge. 

2b) c — i/i, the solution of the Burgers equation has a second point of gradient catas- 
trophe {xi,ti) in correspondence of the trailing edge, there exists a point of gradient 
catastrophe on the M2-branch of the one-phase solution for t > ti and there is a double 
leading edge. 




Xi X ^1 



2a) 2 b) 
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3a) ui < c < 1/2] there are points of gradient catastrophe in the one-phase solution on 

the Mi-branch and M2-branch for t > and there is double leading edge. 

3b) z/3 < c < z/4 ; there are points of gradient catastrophe in the the one-phase solution 

on the Ms-branch and one on the M2-branch and there is a double trailing edge. 




0^0 

3a) 3b) 



4) < c < vs; there is only the point {x — 0,t — 0,u — 0) of gradient catastrophe of 
the solution of the Burgers equation. The curves x~{t) and x'^{t) are C^-smooth for 
all t > 0. 





t 




g=i /-^^ 


x(t)A 








/ 



X 

For c = 1^2, the curve x~{t) looses the C^-smoothness for a single value of t > 0. For 
C — U2 the curve x'^{t) looses the C^-smoothness for a single value oit>0. 

5) < c < i^i; there is a second breakpoint in the solution of the Burgers equation 

for a;i < and ti > 0; there is a point of gradient catastrophe in the one-phase solution 
on the M2-branch for t > ti, there are a double leading edge, a leading-trailing edge 
and a double trailing edge. 
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4.6 Conclusion 

In this paper we have studied the bifurcation diagram of the solution of the Whitham 
equations for a one-parameter family of initial data. On the x — t plane we have 
characterized the bifurcation diagrams in terms of particular singular points which we 
have called double leading edge, double trailing edge, leading-trailing edge and points 
of gradient catastrophe. The x — t plane is split in regions where the solution of the 
Whitham equations has genus g = 0, 1, 2. The genus g = 2 solution survives only for 
a finite time. This result is in agreement with the result in |2^,[^| where it is shown 



that the solution of the Whitham equations with monotone polynomial initial data 
has a universal one-phase self-similar asymptotics. This implies that the large time 
behaviour of the bifurcation diagram of the solution of the Whitham equations can be 



described analitically. Indeed, as follows from p2[, for all times t > T, where T > is 
a sufficiently large time, the solution of the Whitham equations for the one parameter 
family of initial data (|1.8|) is of genus one inside the interval x~{t) < x < x^{t), where 
x^{t) ~ —16.82^5 and x~^{t) ~ 1.58t5. It is of genus zero outside this interval. 



Acknowledgments. I am indebted to Professor Boris Dubrovin who posed me 
the problem of this work and gave me many suggestions for its solution. 
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Figure 4.4: On the plots a) and b) the time of gradient catastrophe of the one-phase 
solution is plotted as a function of c. The circles represent the points of gradient 
catastrophe on the U2 branch, the triangles are the points of gradient catastrophe on 
the M3 branch and the squares are the points of gradient catastrophe on the Ui branch. 
The solid lines represent the time of gradient catastrophe of the zero-phase solution as 
a function of c. 
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